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Analysis of a Class of Likelihood Based Continuous Time Stochastic 
Volatility Models including Ornstein-Uhlenbeck Models in Financial 

Economics. 
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In a series of recent papers Barndorff-Nielsen and Shephard introduce an attractive class of con¬ 
tinuous time stochastic volatility models for financial assets where the volatility processes are functions 
of positive Ornstein-Uhlenbeck(OU) processes. This models are known to be substantially more flexible 
than Gaussian based models. One current problem of this approach is the unavailability of a tractable 
exact analysis of likelihood based stochastic volatility models for the returns of log prices of stocks. With 
this point in mind, the likelihood models of Barndorff-Nielsen and Shephard are viewed as members of 
a much larger class of models. That is likelihoods based on n conditionally independent Normal random 
variables whose mean and variance are representable as linear functionals of a common unobserved Pois¬ 
son random measure. The analysis of these models is facilitated by applying the methods in James (2005, 
2002), in particular an Esscher type transform of Poisson random measures; in conjunction with a special 
case of the Weber-Sonine formula. It is shown that the marginal likelihood may be expressed in terms of 
a multidimensional Fourier-cosine transform. This yields tractable forms of the likelihood and also allows 
a full Bayesian posterior analysis of the integrated volatility process. A general formula for the posterior 
density of the log price given the observed data is derived, which could potentially have applications to 
option pricing. We also identify tractable subclasses, where inference can be based on a finite number of 
independent random variables. We close by obtaining explicit expressions for likelihoods incorporating 
leverage. It is shown that inference does not necessarily require simulation of random measures. Rather, 
classical numerical integration can be used in the most general cases. 


1 Introduction 

Barndorff-Nielsen and Shephard (2001a, b)(BNS) introduce a class of continuous time stochastic 
volatility (SV) models that allows for more flexibility over Gaussian based models such as the Black- 
Scholes model[see Black and Sclioles (1973) and Merton (1973)]. Their proposed SV model is based 
on the following differential equation, 

(1) dx*(t ) = (g + (3v(t))dt + v 1 / 2 (t)dw(t) 

where x*(t) denotes the log-price level, w(t) is Brownian motion, and independent of w(t), v(t) is a 
stationary Ornstein-Uhlenbeck (OU) process which models the instantaneous volatility. This model 
is an extension of the Black-Scholes or Samuelson model which arises by replacing v with a fixed 
variance, say cr 2 . The additional innovation in BNS is that modeling volatility as a random process, 
v(t), rather than a random variable, not only allows for heavy-tailed models, but additionally induces 
serial dependence. This serial dependence is used to account for a clustering affect referred to as 
volatility persistence. The work of Carr, Geman, Madan, and Yor (2003) discuss this point further. 
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See also Duan (1995) and Engle (1982) for different approaches to this type of phenomenon. The 
model of BNS has gained a great deal of interest with some related works including Carr, Geman, 
Madan, and Yor (2003), Barndorff-Nielsen and Shephard (2003), Eberlein (2001), Nicolato and 
Venardos (2001), Benth, Karlsen, and Rcikvam (2003). See also the discussion section in Barndorff- 
Nielsen and Shephard (2001a). 

One current drawback of this approach is the unavailability of a tractable analysis of likelihood 
based stochastic volatility models for the returns of log prices of stocks. These models are based on 
the integrated volatility process r(i) = f 0 v(u)du. Several MCMC procedures have been proposed 
to handle subclasses of these models requiring simulation of points from random processes. See for 
instance, Roberts, Papaspiliopoulos and Dellaportas (2004) and the discussion section in Barndorff- 
Nielsen and Shephard (2001a). 

In this paper, we shall actually provide a complete analysis of a significantly more complex 
class of likelihood models. Specifically models where r is expressible as a linear functional of a 
Poisson random measure. This includes the superposition processes mentioned in Barndorff-Nielsen 
and Shephard (2001a, b) and much more general spatial models for r. Our results are therefore 
applicable to a wide range of models and applications. We first present a description of these models 
similar to the framework outlined in Barndorff-Nielsen and Shephard (2001a). We then show how 
the results in James (2005, 2002) are easily applied to this setting via the usage of a Bessel integral 
identity involving the cosine function. That is, a special case of what is called the Weber-Sonine 
formula. This leads to an interesting series of tractable characterizations of such processes, including 
an identification of simple subclasses of these models. As a byproduct, we derive the posterior 
predictive density of what could be considered as models for the log price of stocks. This may 
prove useful to applications in option pricing. Moreover, our methods do not require simulation 
of random measures and in the most general cases can be handled by more classical numerical 
integration methods. We point out also that procedures to simulate from random measures often 
require explicit knowledge of the Levy density associated with an infinitely divisible random variable. 
This is an important point as there are some interesting cases where the probability density of a 
random variable is known explicitly but its corresponding Levy density is unknown. Our results show 
that one only requires knowledge of the form of the Levy exponent or log of the Laplace transform of 
a corresponding random variable. Inference using the classes of models that we identify in sections 
2.5 and 4 can be performed based on at most J < oo independent latent random variables. In the 
case of section 4, J = n + 1, where n denotes the number of observations. Section 2.6 discusses 
another class where exact calculations of a different form are easily obtained. Section 5 describes 
exact expression for likelihoods of generalized types of leverage effects models. 


Remark 1. The appearance of integrals involving Bessel functions is certainly not new to 
applications in finance as can be seen in the case of the important work of Yor (1992) on Asian 
Options. See also Carr and Schroder (2004). We shall however employ a different, but certainly 
related, integral identity. 


1.1 Likelihood model and representation 

The model of Barndorff-Nielsen and Shephard (2001a, section 5.4) translates into a likelihood based 
model as follows. Let JQ for i = 1,... ,n denote a sequence of aggregate returns of the log price 
of a stock observed over intervals of length A > 0. Additionally for each interval [(i — 1)A, zA], let 
Tj = r(?'A) — r((i — 1)A). Now the model in QJ implies that A, jr^, /?, fi are conditionally independent 
with 


( 2 ) 


Xi = llA + Tifi + T 
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where tj are independent standard Normal random variables. Hence if r depends on external pa¬ 
rameters 9, one is interested in estimating (//, (3,9) based on the likelihood 


( 3 ) 


-Sf(X| y,{3,9) 



Y[(j}(Xi\lJ,A + pTi,Ti) 
_i= 1 


/(t 1,... ,T„|0)dTi, ...,d,T n 


where, setting Ai = (X,; — yA), and A = n 1 

+ 0Ti,Ti) = e Ai P-^Tr 1/2 e- A W^e- T *P 2 / 2 

\/2ir 

denotes a Normal density. 

The quantity /(n,..., r n |0) denotes the joint density of the integrated volatility based on the 
intervals [(* — 1)A,*A] for i = 1 Barndorff-Nielsen and Shephard (2001a) note that the 

likelihood is intractable and hence makes exact inference difficult. The apparent intractability is 
attributed to the complex nature of /(ti ,..., r„ |0) which is derived from a random measure. However, 
we shall show that in fact it is quite easy to deal with /(ti, ..., t„| 0) for more general r by means 
of the Poisson partition calculus methods outlined in James (2005, 2002). Rather, the stumbling 
block which currently prevents one from integrating out the infinite-dimensional components in 
the likelihood, is inherent from the Normal distribution of A ?: |r,;, /3, //. Quite simply the Normal 
assumption yields exponential terms of the form 

Q- A V( 2 Ti) ra p ier than e~ TiAi . 


In the next section we apply an integral identity to circumvent this problem. 


1.2 Bessel integral representation of the likelihood 


In order to calculate © we first employ a Bessel integral identity which we state in more general 
terms. Suppose that J v (x) denotes a Bessel function of the first kind of order v. Then for v > —1, 
and numbers a,p 

pOO 

(4) p~ 2(«+i) e -« 2 /V = 2 v+1 a ~ v / J v {at)t v+1 e~ pH2 dt. 

Jo 

This is a special case of the Weber-Sonine formula. See for instance Andrews, Askey and Roy (1999, 
p.222) and Watson (1966, p. 394 eq. (4)) for the identity and also those references for Bessel functions. 
Taking for each i, p 2 = t^/2, a = |Aj| and v = —1/2, it follows the marginal likelihood is given by, 


( 5 ) 


where 


if(X|/x,/3,0) 


a nAf) 


E 


IT 


,-(2/-/2+/3 2 /2)n 


n 

Y[cos(yi\Ai\)dyi 

i=1 


E 


n« 

i=i 


-(l/?/2+0 2 /2)T 4 


n n 

/ U^ / 2 + 02 / 2 ) ^f^ ■ ■ • > T «I 6 ') II dn. 

J n i =i 1=1 


Remark 2. Notice that the expression © has nothing to do with the distributional properties 
of r. The appearance of the cosine in © is due to the identity 


(- 1 )' 


X 


^ o T(k + l/2)k\2 2k ~ 1 / 2 


= X 1 ! 2 J_ 1 / 2 ( 2 ;) = \ — cos(x) 


where J_ i/ 2 (*) is a Bessel function of the first kind of order —1/2. See for example Andrews, Askey 
and Roy (1999, p. 202). In other words we are using the classical Fourier-Cosine Integral 


1 

7 r 


cos(y|Aj|)e 2 ‘ dy = 
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2 Evaluation of the likelihood for general r 

The representations in © and allows us to immediately apply the results in James (2005) to 
obtain a full analysis for quite general r, which we now describe. Let N denote a Poisson random 
measure on some Polish space Y with mean intensity, 

E[N(dx)\v} = v(dx). 

We denote the Poisson law of N with intensity v as P(dN\v). The Laplace functional for N is defined 
as 

E[e -JV( ^' ) |i'] = f e~ N ^^P(dN\u) = e _A ^ 

Ju 

where for any positive /, N(f) = Jy f(x)N(dx) and A(/) = / r (l — e~^ x>> )u(dx). M denotes 
the space of boundedly finite measures on Y [see Daley and Vere-Jones (1988)]. We suppose that 
Ti = for i = 1,..., n where are positive measureable functions on Y. Notice now 

that the index i = 1,... ,n need not correspond to fixed intervals involving A. With this in mind, 
let (wi ,..., w n ) denote arbitrary non-negative numbers. Define for i = 1,... ,n, functions Ri{x) = 
Ej=i w jfj ( x ) an< l VRi{dx) = e~ Ri( ^ x ' > v{dx). Then all our results will follow from the following special 
case of James (2005, Proposition 2.1), which can be viewed as an Esscher-type transform, 

( 6 ) e ~ N G2i=i Wif df>(dN\v) =F(dN\v Rn )e- A ^=i Wifi) . 


Additionally the following decomposition is sometimes useful 


E 


3 -^(Er=i | 


= =E 


^-N(w i/i) I 


n 

n 4 


,-N(wifi) 


VRi_ i 


i=2 


This expression appears in James (2002) and may be obtained by repeated application of ©. 

Now, throughout, for each n > 1, define D n (x) = E™=i( 2/i/2 + /3 2 /2)fi(x). This is a special case 
of Er=i Wi fi with Wi = (yf/2 + (3 2 /2) for i = 1,..., n. The following result is immediate from an 
application of Fubini’s theorem, m and the representation ®. 

Theorem 2.1 Suppose that Ti = N(fi) for i = 1,... ,n where N is a Poisson random, measure on 
Y with intensity v. Then setting Wi = yf /2 + (3 2 /2 for i = 1,..., n, in ©, the likelihood m can be 
expressed as 

p nA 0 t- n 

Jz?(X|/x,/?,6») = —— / g- A (f2n) TTcos(yi|Ai|)ch/i. 

17 M i= l 

□ 


2.1 Posterior distribution of parameters 

Theorem 2.1 shows that Bayesian inference for (fj.,/3,6) may be described as follows. 


Proposition 2.1 Suppose that t depends on a d-dimensional parameter 9. Then if q(d9), q(df3), 
q(dg) denote independent prior distributions for (/?, p, 9), their posterior distribution can be written 


q{dp,dy,de |X) oc 


-[A„(Q„)-? 


At3] q{d9)q{dfd)}Y[cos{y l \A l \)dyiq(dy), 


where Ag denotes the dependence of A on #.□ 
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2.2 Posterior distribution of the process 

The above results describe the behaviour of the finite-dimensional likelihood and parameters. It is 
useful to also obtain a description of the underlying random process given the data. This allows 
one to see directly how the data affects the overall process. Moreover, combined with the results in 
James (2005), it provides a calculus for more general functionals. Define the measure, 

n 

^x(dy) = ^ e -^n)-nA 0 ] JJ c os( y< |^|)dy i /JSf’(X| A i, J 9,6l). 

For notational simplicity we suppose that (/j, (3, 9) are fixed. The next result also follows immediately 
from an application of Fubini’s theorem, © and the representation ©. 

Theorem 2.2 Suppose that the distribution of X is given by w, and that r and N are defined 
by the specifications in Theorem 2.1. Let D. n (x) = JT_i(?/f/2 + (3 2 /2)fi{x). Then the posterior 
distribution of N |X is given by the mixture 

[ P(dN\vn n )£x(dy) 

J 

which determines the posterior distribution of t and related quantities. P(ci/V|t'n n ) can be viewed 
as the posterior distribution of N given the information in Y, X and corresponds to the law of a 
Poisson random measures with mean intensity 

va n (dx) := e~ n ^ x) v{dx) = iy(dx)e~^^ y ^ 2+ ^/ 2 '> fi ^. 


□ 


The next result, which gives an expression for the posterior Laplace functional of N, is an immediate 
consequence of © combined with Theorem 2.2. 

Proposition 2.2 The posterior Laplace functional of iV|X, according to the Theorem 2.2, is given 
by 


e - NU) P{dN\vn n 

M 

for f such that A(/ + fi n ) < oo. □ 


^x(dy) = [ e-^ +n ^- A ^J2x(dy) 

JR 1 } 


2.3 A general posterior predictive density for the log price 

We now define a random variable similar to © which can be thought of as representing the log-price 
and give an explicit expression for its posterior density given X. The random variable is defined as, 

(7) X = fiK + f/3 + f 1 / 2 e 

where A denotes a general positive quantity, e is a standard Normal random variable independent 
of all other variables, and f = N(f) for some positive function / such that Laplace transform of f 
exists. 

Proposition 2.3 Suppose N and the data structure of X is defined as in Theorem 2.2. Let X be 
defined by 0- Denote its marginal density as f-y(-\P,/T) and its posterior density given the data X 
from m as fy (’|/3, X). Then the following results hold 

(i) fx(x\/3,p) = I e (-/^)/3 J- e- A ^/ 2+ ^^^co S (y\x-pA\)dy. 
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(ii) The posterior density of the log stock price given X is, f^{x\j3, p, X), given by 


7T 




e -[A(u) n+1 /+fi»0-A(fi„)] cos (y\ x 


liA\)dy 


^x(dy), 


where w n+ \ = (y 2 /2 + p 2 /2).0 


Proof. Setting ri ra+ i(a;) = Ll n (x) + w n+ \f{x), the results follow from © and Theorem 2.2, 
using the fact from © that, 

e -^(™ B +i / ~)p(diV|r'fiJ = P(d7V|z/ f2ri+i )e-[ A ^”+ l/+n ’* ) - A(n " ) ]. 


□ 


2.4 Simplifications for a class of r via an inversion formula 

We have shown that for r modeled quite generally that its contribution to the likelihood m is only 
through the exponent AQT)" =1 Wiff). That is through the form of YAi=i w ifi an d v - With a view 
towards choosing r which are the most tractable we present the following interesting result. 


Theorem 2.3 Suppose that for arbitrary non-negative (wi ,... ,w n ), and an integer J, there is an 
array of non-negative numbei's ( aij ) such that A(^" =1 Wifi) = % A(E"=i WiCbij]hj). Where (hj) 

are non-negative functions on "V such that A(u>hj) < oo for all u > 0. Let {T\,... ,Tj) denote 
J independent random variables with respective Laplace transforms E[e~ ujTj ] = e~ A ^ ujhj ' > for j = 
1,..., J. Moreover „S?(X| p, (3,9) denotes the likelihood given in ©. Then, 


(i) E [e-E^m] = e- A (Er=i^“oi^) = J] j =1 E[e 


r-EI 


i]Tj 


(ii) Jgf(X|/x, 0 , 0 ) = E [nili + 0\£ J j= 1 E/=1 OijTj ]) 

respect to the distribution of (Ti,..., Tj). 


, where the expectation is 


□ 


PROOF. Statement (i) is immediate from the specification of AQT)” =1 Wifi). Statement (i) implies 

that that one may replace YAi=i w i T i with Y17= i ®«E/=i a ijTj]- Now setting each wi = (y 2 / 2+ f3 2 /2) 
for i = 1,..., n, one uses © and © to conclude the result. □ 

Statement (ii) of Theorem 2.3 allows one to approximate the likelihood by the simulation of J 
independent random variables. It also demonstrates that it is rather straightforward to conduct para¬ 
metric Bayesian or frequentist estimation procedures, where (Ti,..., Tj) are viewed as independent 
latent variables. The next proposition puts this in a Bayesian framework. 

Proposition 2.4 Suppose that (Ti,...,Tj) depend on external parameters, say 9. Then assuming 
a joint prior q(d9,dp,dP), posterior inference may be obtained based on the model derived from 
augmenting the likelihood in Theorem 2.3. That is, the joint distribution of (X, T \,..., Tj, 9, /a, /3) 
given by 

n J J 

II 'HA-, //A + 0[J2 U,E a d T j}) 

i=1 j =1 3=1 


Ufr 3 (T,) 


3=1 


q(d9 , dp , d0). 


□ 
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3 Tractable expressions 

It is noted that, at first glance, one may find it difficult to work with the expressions involving cosines. 
Here, influenced by some arguments in Devroye (1986a), we give a representation of the likelihood 
that can be numerically evaluated via the simulation of random variables. First let p = {p \,... ,p n } 
denote a vector of positive numbers and for each *,let 


H(Vi\Pi) = 


y/2irp~i 


e 2p i for yi > 0 


denote a half normal density. Now, notice that 0 < 1 — lT=i C0S (Vi) < 2, and 


( 8 ) 


1 - H c °s(2/ii^i 


2=1 


Ef=i A fpi „ 

H(yi\pi)dyi = 1 - e 2 = C n { A, p) 


From these facts we describe a joint density 

Proposition 3.1 Augmenting the expression in 0) leads to a joint density of an array of positive 
random variables Y = {Yi n ,..., Y n , n } given by, 

r fix [! - IT=1 COS^ilAil)] 

n(y|Pj_ C„(A, p) 

Equivalently, for k = 1 the conditional density of Yh in \Yi >n , ■ • •,bc-i,n is proportional to 

yy n Pi 1 

[1 - \ k cos(y k \A k \)\H(y k \p k ), where X k = e ~Yi=k+i 1L=1 cos(yi|A|) for k = 2, 
e~Yi =2 —, and \ n = nEi cos(j/j|A|). 

Define the function 


,n- 1, Ai = 


T„(/?,0) := 


1 


=,-A(n„) 


71 JR* 


JJ dyi = E 


i=1 


,~0 2 Ti 




< E 


n 

.2 — 1 


These points lead to following representation of the likelihood. 

Proposition 3.2 Suppose that for fixed n, E JT = i j/F < °°’ then the likelihood in Theorem 2.1 
may be written as 

Af3 r rDA ni r P -H^n 


r n (p,e)- Cn(A J p) E 


nr =1 


where the random vector {Yi Y n „} has its joint distribution described by proposition 3.1, and 

nn{x)=Y.U[{Y?, n + P 2 )/i9i{x) 

Remark 3. Proposition 3.2 shows that one may approximate the likelihood by simulating 
random variables described in Proposition 3.1. Such an approach should work well with a Bayesian 
procedure for estimating the parameters (y, (3, 9). Methods to easily sample the random variables in 
proposition 3.1, may be deduced from Devroye (1986a, b). Alternatively one may sample from the 
densities H(yi\pi). One may also use other densities. 


4 Analysis of the BNS-OU model 

In this section we will show how our results apply to the basic integrated volatility model of Barndorff- 
Nielsen and Shephard (2001a, b). We shall refer to this model as the BNS-OU model. First suppose 
that A is a Poisson random measure on (0,oo) x (—cx),oo) with intensity 


u(du , dy) = p{du)dy 
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where p is the Levy density of an infinite-divisible random variable, say T, with Laplace transform 
for to > 0, 

r-OO 

E[e"“ T ] = e-'/’M where ip(co) = (1 - e~“ u )p(du). 

Jo 

Now we model the background driving Levy process (BDLP), say z, as a completely random mea¬ 
sure which is expressible in distribution as z(dt) = f 0 uN(du,dt). Note that for any non-negative 
function g on (— 00 , 00 ), it follows that 

/ OO 

g{y)z(dy ) = N{f g ) 

-00 

where f g (u,y) = ug{y ) on (0, 00 ) x (— 00 , 00 ). Additionally, 


E 



e~ N{f ^P(dN\iy) 


/M 


_ e - f- 00 i’( 9 (y))dy _ e -A(/„). 


One may express the Barndorff-Nielsen and Shephard (2001 a, b) integrated OU process r as 


(9) r(t) = A ^(l-e At ) [ e y z(dy) + f (1 - e A(t v) )z{dy)} 

J —00 J 0 

where n(0) := Vo = e y z(dy). The form in 0 is taken from Carr, Geman, Madan and Yor (2003, 
p. 365). It follows that for any s <t, [r(t) - r(s)] = z(g Sit ) = N(f s>t ) where f s , t {u,y) = ug Sjt (y) and 
A g s ,t(y ) equals, 


(10) e" As (l - e - A (‘- s V^<o } + (1 - e- x ^)I{ s<y < t} + e ~ Xs (l - e - A (‘-«)) e A vj {0<1 ,< B} . 

The first component in m represents the contribution from vq. Specializing this to s = (i - 1)A 
and t = iA one has n = z(g itl + g it 2 ) = N(fi) where ffiu, y) = u[g itl (y) + ft, 2 ( 2 /)] and 


(11) g iA (y) = A" 1 [(l - e- A ( iA -^)/ {(i _ 1)A<3/ < iA} + - e - AA )e^ { ,< 0} ] 


and 


(12) ft, 2(y) = A _1 e _A ^ _1 ^ A (l - e _AA )e Ay / { o< !/ <(j_i) A} . 

Now for i = 1,..., n, set r* = A _1 [J]^ =i 'U)fce _A ^ fe_1 ^ A ](l — e _AA ). Now notice that for any sequence 
of numbers, the simplest expression will be obtained by utilizing the following facts. 

n 

(is) ^ 2 ^j\gj,i(y) + 93Av)} = r ^ v for y < o 

j =1 

and for i = 1,..., n 

n 

(14) ^2wj[g jtl (y) + gj,2(y)] = C(j/K, A+i) for (i - 1 )A < y < iA. 

i=i 

Where for each i, ((y\wi,ri + i) = [A -1 rft(l — e - A 0 A -y)) + ri+ie Ay ]. 

Proposition 4.1 For 0 < s < t, let r(t) — r(s) be defined by 0) and m Then the results of 
Theorem 2.1 and 2.2 hold with fi(u,y) = u[gi,i(y) + ft, 2 ( 1 /)], ift = {yf /2 + (3 2 / 2), as described 
in m and m- In particular, using a change of variable, 
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(i) e~ A ^*= lWi ^ = e - ^ 0 ^ 7 " 1 ^e~^ n ^ Wn ^ [j"-/ e~^ i ^ Wi ^ ri + 7> 

(ii) $(wi\r i+1 ) = fi_ e -x a X~ 1 ip(r i+ ie XzA (l - u) + X~ 1 w i u)^, for i = 1,..., n - 1 
(in) $(w„) = e _ AA 

(iv) $ 0 (^ 1 ) = fg V’( r i u )if; w/iere e _<I>( A ri ) = E[e _rit ’ 0 ]. 

□ 


Remark 4. Expressions of the form in [(iii)] of Proposition 3.1 are known to be a key component 
in option pricing using the BNS-OU model. However explicit calculations have only been given for 
a few cases. See Barndorff-Nielsen and Shephard (2003), Nicolato and Venardos (2003) and Carr, 
Geman, Madan and Yor (2003). Note that if for y > 0, we change the Lebesque measure, dy, to 
e Xy dy , the calculations for $(w*|ri+i) for i = 1 ,...,n, where $('ui rl |r rl+ i) = 4>(ui n ), are greatly 
simplified. See section 4.2 for a closely related discussion. 


5 Analysis of a simple class of models 


This last section, which is based on a class of models from section 2.4, examines models which are 
the most tractable and we believe still flexible enough to be applied to general classes of problems. 
Implicitly, we are taking a Bayesian nonparametrics viewpoint of seeking random measures as priors 
which are both flexible in a modeling sense and easily manipulated. For concreteness, we start out 
with a variation of the Barndorff-Nielsen and Shephard (2001a,b, 2003) integrated OU process r. 
Here we set, 

(15) r(t) = A _1 [(l — e _At ) f e v z(dy)+f (1 - e~ v )z(dy)] : 

J — 00 J 0 

where again t>(0) = Vo = e v z{dy). Interestingly from Barndorff-Nielsen and Shephard (2003, p. 
282), one has the following distributional equivalence of marginal distributions, 

pXt /*A t 

A -1 / (1 — e~ v )z(dy ) = A” 1 / (1 — e~ x< ' t ~ v ^)z(d\y) 

Jo Jo 

where the right hand side equates with the model in section 3. However, now for s < t, 

A[r(f) - t(s)] = e _As (l - e _A(t_s) ) [ e v z(dy)+ [ (1 - e~ y )I {Xs<y < xt} z(dy). 

J —00 J 0 


Specializing this to t = iA and s = (i — 1)A yields 


Ti = 



9 i{y)z{dy) 


z{9i) = N{fi) 


where fi(u,y) = ug.Jy) with 

9i(y ) = A -1 [(l - e- v )I {X (i-i)A <y <\iA} + e" A(l_1)A (l - e~ XA )e y I {y < 0} ]. 


Now notice that for any sequence of numbers, the simplest expression will be obtained by utilizing 
the following facts. 


J2 w j9j(v) = A-^^e-^-^Kl - e _AA )e y for y < 0 

j=1 *=1 


( 16 ) 
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and for i = 1,..., n 
(17) 


^ ' / Wjgj(y ) = A 1 Wi( 1 — e v ) for X(i — 1)A < y < XiA. 
l=i 


More generally suppose that for t > s, r(t) - r(s) = z{g Stt ) = N(f Syt ) where f Stt (u,y) = ug S:t (y) 
and 


(18) 


g s ,t(y) = h ltSt t,x(y)I{x s<y <xt} + h 2 , s ,t,xF(y)I {y < 0} , 


for hi s t} x(y) and F(y) non-negative functions satisfying suitable integrability conditions and h 2 , s ,t,X 
a positive quantity not depending on y. Hence, one could choose for each i, 


(19) 


5^(l/) — h , l,X,i{jj')I{X(i— l)A<y<AiA} T ^2,A,iT\y)-Z{y<o} > 


for arbitrary positive functions h\ } x,i, F and h 2 ,x,i whose form is determined by the general difference 
r(t) — r(s) for t > s. These models all exhibit behavior similar to 111 (ill and O- That is for any 
sequence of numbers (w ±,..., w n ), it follows that still Y^i=\ w i T i = -KSiLi w i9i) = -W'QZiLi 
for fi(u, y) = ugi{y) and gi now given by 11811 . Additionally, the most important feature is preserved. 
That is, 


( 20 ) 


= \^Z w i h 2 ,x,i]F{y) for y < 0 
i=i 




and for i = 1,... ,n, J2]=i w j9j(y) = Wihi,x,i(y) for A (i - 1)A <y < XiA. 
Now let N denote a Poisson random measure on (0, oo) x (— 00 , 00 ) with 


( 21 ) 


u(dw, dy) = p 1 (dw)i] 1 (dy)I{ v>0 } + p 2 (dw)p 2 (dy)I{ y < 0 }, 


where p\ and p 2 are Levy densities generating Levy exponents ifi and ip 2 , and pi, p 2 are non-negative 
sigma-finite measures. It follows from that 


( 22 ) 

where 


il’i(^2' w i9i{y))Vi(dy) + / ip 2 (^2wigi(y))p 2 (dy) = $ 0 (s„) + ^$i(wi) 


2=1 


2=1 


2=1 


piA nO n 

$i(wi) = / il’i(w i hi,i i \(y))r]i(dy) and $ 0 (s„) = / Wjh 2 ^x\F{y))p 2 {dy) 

J \(i—l)A J — 00 „_i 


for i = 1,..., n, where s n = X)"=i w i^,i,X- In the case of ll(ill and for pi = p 2 , and pi(dy) = 
p 2 (dy) = dy, s n = E"=i 'U'ie _A(l_1)A ](l - e _AA ). One has for * = 1,..., n, 


(23) 


(i_ e ,iA 

$i(wi)= / i/j(X~ 1 w i u )——= / ^{wiX -1 {l - e~ v ))dy 

J{ i_e-A(i-i)A) i-w Jx(i-i)A 


XiA 
A(2-1)A 


and *I>o(s rl ) = f* ip(X 1 s n u)- is the Levy exponent corresponding to the prior distribution of vo 
evaluated at s n . 

Theorem 5.1 Let N denote a Poisson random measure on (0, 00 ) x (— 00 , 00 ) with intensity v 
defined in m- Define r by HE) and UB and the general specification of v above. Let <$>j for 
j = 0,1,..., n denote the quantities defined by 1221 ) . Let (vo, T \,..., T n ) denote n+ 1 independent 
random variables with respective Laplace transforms for i = 0,1,. .., n. Then, 



Lancelot F. James 


11 


(i) E[e-£"=i“ iT *] = nr=i =E[e- s ”^°)|^nr=i E [ e_u ' iTi l i/ ] 

(ii) Jf(X\p,P,9) = [Yli =1 f™ (/>{Xi\fj,A +/3zi,Zi)f Ti {zi-biy)dzi] f Vo (y)dy where bi = h 2 ,i,x 
fori = 1 ,,n. /t i; • ■ •, /t„) denotes the densities of the corresponding random variables .□ 


Proof. 

3i].n 


The result is a special case of Theorem 2.3 where one replaces w i T i with X]”=i w i [biVo+ 


Bayesian inference may be conducted using the following result. 


Proposition 5.1 Suppose that (T \,... ,T n ,v o) depend on external parameters, say 9. Then assum¬ 
ing independent priors q(d9), q(dp) and q(d/3), posterior inference may be obtained based on the 
model derived from augmenting the likelihood in Theorem 5.1. That is, the joint distribution of 
(X, Ti,..., T n , v 0 , 9, p, /?) given by 


□ 


Yl^XilnA + liZ^Zi) 

_i= 1 


q(dp)q(dp) 


n fM 

_i =1 


biV o) 


fv 0 (vo)q(d9). 


5.1 Predictive density of the stock price 

The model m suggests for any time interval [s, t] that the return of the log stock price, say X s p, is 
given by the model, 

(24) X 8 , t = p(t - s) + [r(t) - t(s)]/3 + [r(t) - T(s)] 1/2 e s , t 


where e Sy t is an independent standard Normal distribution. Our results yield an explicit tractable 
expression of a predictive density of X Sit given previously observed data X 


Proposition 5.2 For t > s > nA, let X Sjt be defined according to \%4\) - Let t be defined by the 
general specifications in Theorem 5.1 Let fx ± ,..., /r„, fv 0 denote the densities of the corresponding 
independent random variables. Then the predictive density of X s ^\X\,... ,X n is given by the formula, 



(f(x\p(t 


s) + /3z Stt ,z S! t)fT a , t (z s ,t - b Stt v)dz Stt 


r(X\v)f Vo {v)dv 


where, r{X\v) = [JlLi fo° 0(X,:|/iA + (3 z t , Zi)f Ti {zi - biv)dzi] /Jgf(X| p,f3,9) andb Sft = h 2 , s ,t,\- The 
quantity, fr, t denotes the density of an independent random variable T S}t with law determined by 

the Levy exponent $ s ,t(w) = J x * ipi(whi ti ,\(y))rii(dy). □ 


PROOF. The result follows by noting that )C”=i w i9i{v) + w n+igs,t(y ) has the same structural 
form as 12(J|| . Where g S:t is defined in (HU . □ 

5.2 Example 

The expressions in Theorem 4.1 suggests that an easily analyzed model would arise if (vq, T\ ..., T n ) 
were all from GIG class of densities. Here, going back to the variant of the Barndorff-Nielsen and 
Shephard model characterized by we shall show that the choice of a stable law yields very nice 
results. Recall that the Levy exponent of a stable law of index 0 < a < 1, is such that = u> a /a. 
Recall also that the case of a = 1/2 leads to the Inverse Gamma distribution of index 1/2, and that 
the Inverse Gaussian arises from an exponential tilting of this law. Using this fact we arrive at the 
following result. 
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Proposition 5.3 Suppose that r is specified by 1751) with v(ds,dy) = s -a-1 /[r(l — a)]dsdy for 
0 < a < 1. Then the random variables (vo, T) ..., T n ), appearing in Theorem 5.1, are independent 
stable random variables of index a. The respective Levy exponents are, 

(i) $oM = w Q A-“(l - e~ XA ) a /a and 

(ii) $i(w) = io a \~ a fiZe-Hi-i) a fori = l,...,n 

(Hi) If one instead uses r]i(dy) = e~ v dy, then 

= - T ^ ZTT - A " a [( 1 - e - AlA ) Q+1 - (1 - e _A ( i_1 ) A ) a+1 ] 

a(a + 1 ) 

for i = 1,. ..,n. 

□ 


Notice that the proposition above shows that a change from Lebesque measure to rji(dy) only 
affects the constants in the Laplace transform and generally preserves the distributional property 
of the (Ti). This fits into what has been evidenced in Bayesian nonparametric problems where the 
choice of quantities such as rg, 1)2 are done mainly for computational convenience. The rationale is 
that viewing the specifications for r as a prior model, experience from the Bayesian nonparametric 
literature suggests that many such choices of r will eventually lead to the similar conclusions in the 
presence of enough data X. Note however that we do not advocate removing the dependence of the 
Levy exponents on (i, A), as this is related to the data. We close by noting it is always possible 
to arrange for the random variables (vg, Ti,..., T n ) to be self-decomposable by choosing if\ and 
if 2 , rji and r ]2 such that the random variables are Generalized Gamma Convolutions (GGC). See 
Thorin (1977) and Bondesson (1979, 1992) for this rich class of models. That is $i(w) = ln(l + 
w/y)ty(dy), for aj,Cj depending on A. %, with ^(0) = 0, is called a Thorin function or measure. 
This applies more generally to the models described in section 2.4. We shall leave it to the reader to 
investigate which of the general models discussed in section 2 are most suitable to their particular 
application. 


6 Extension: SV Likelihood models with correlated jumps in 
price, leverage effects models. 

Recall that the Barndorff-Nielsen and Shephard(2001a,b) OU process v(t), which models the instan¬ 
taneous volatility, satisfies the differential equation 

dv(t) = —A v(t) +dz(Xt), 

where the process z is defined in section 3, and hence the volatility possesses jumps. An important 
extension of the model in m and hence to our general framework described in section 2, is where 
one includes jumps in the log-price model which are correlated with the the volatility v. These types 
of continuous time models, wherein Duffle, Pan and Singleton (2000) is an early reference, serve to 
incorporate the leverage effect discussed in for instance Black (1976) and Nelson (1991). We shall 
be rather brief on this growing literature and refer the reader to the works of Eraker, Johannes 
and Poison (2003), Duan, Ritchken and Sun (2004) and Duffle, Singleton and Pan (2000), for more 
extensive background and rationale for these types of models and its parametric variations. 
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6.1 BNS-OU SV model for leverage effects 

Barndorff-Nielsen and Shephard (2001a, eq. 8) describe this type of extension as follows, 


(25) dx*(t ) = (p + / 3v(t))dt + v 1 / 2 (t)dw(t) + p(d[z(Xt) — E[d.z(A(i)]) 

assuming of course that E[z(At)] < oo. One can incorporate modifications to relax this condition. It 
follows that obviously the log price and the volatility are negatively correlated if p < 0. Thus modeling 
the leverage effect that a fall in price results in an increase in future volatility. Barndorff-Nielsen 
and Shephard (2001a, section 4) discuss further details of this model. The likelihood model based 
on was not explicitly discussed in that paper, and moreover is considered even more challenging. 
However as we shall show, this extension and a variety of natural extensions of the models described 
in section 2 , incorporating a leverage type effect, are easily handled by the type of methodology we 
have presented so far. 

First, assuming a similar framework as in section 1.1, and using the BNS-OU model described 
in section 3, note that 

cOO nOO 

Zi ~ z(g i}3 ) = z(iA\) - z(X(i - 1)A) = / / I{(i-i)A\< y <\iA}uN(du, dy), 

Jo Jo 

where N is a Poisson random measure with intensity v{du 1 dy) = p(du)dy on (0, oo) x (—oo, oo), and 
9i,s{y) = I{(i-i)A\<y<\iA}- Assuming a finite first moment, one has E[zj] = A J Q up(du). Hence 
the model (|25|l implies that A, : |t;, Zi, p. p are conditionally independent with 

(26) Xi = pA + nP + r^ 2 ei + p(zi - E[zj]), 

1 /2 

which may be rewritten for each i , as Xi = (pA — pE[ 2 j]) + (jiP + pZi) +r t e*. Hence on may write 
the expression in as, 

(27) Xi = (pA - p&[zi\) + (z{gi,i + gi^)P + pz(g it 3 )) + \jz{g,p + gi, 2 )^i, 
which obviously may be further expressed in terms of a common Poisson random measure. 


6.2 A General class of likelihoods which incorporate leverage type effects 

We note that from our point of view a model such as m poses no additional complications. Similar 
to section 2 , we will obtain exact expressions for likelihoods of quite general extensions of models 
with correlated jumps in price and volatility. As before, for a general Poisson random measure on 
with intensity z/, let r* = N(fi) for i = 1,... ,n. Additionally, for real-valued functions ipi,... , ip n on 
Y, each satisfying the condition A(|^i|) < 00 , for i = 1,..., n, define, 7 ; := N(ifi). Now a general 
version of (E3) is given by the case of conditionally independent 

(28) Xi = pA + p^p + nP + T^ 2 £i. 

Let y = R x R + , and assume that N depends on a parameter 9. Then the likelihood of X| p, /?, 9 , p 
determined by 1251) can be expressed as 


(29) if(X| p,p,9,p)= f 


4>{Xi\pA + P7i + pTi , Ti) 


.»= 1 


/((n, 71). • • •, (Tn, ln)\0) d'yidTi 


i= 1 


where, 


(j)(Xi\pA + Pll + p Ti , Ti) = 1 r -V 2 [ e -(^-™) 2 / (2r i ) ]e -r^/2 

V 27r 
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The difficulty in evaluating this likelihood now manifests itself in the term in brackets where both 
7 j and r, are functionals of a common Poisson random measure, and moreover are not pairwise 
independent across i. Clearly one could apply 0 however the cosine representation, now involving 
random terms, would generally lead to expressions which are less aesthetically pleasing. Here we will 
simply use an identity deduced from the characteristic function of a Normal random variable. This 
is equivalent to 0 - We close by describing the explicit likelihood. 

Theorem 6.1 Suppose that N is a Poisson random measure with intensity v on 'V. Then define 
the complex valued function T n (x) = ^"=1 [p(/3 + £yi)]</ 7 :(x) on Y, where £ denotes the imaginary 
number. Then defining fl n (x) as in Theorem 2.1 with Wi = yf/2 + /3 2 /2 for i = 1 ,...,n, the 
likelihood [Mi can be expressed as 

p nA/3 r n 

^(X| M , /J, e , p) = —/ e- A ( n " +T ") n S A ^d Vi . 

(^ 7r ) JR” i=1 

where (ipf) and ( fi) are chosen such that A(f2 n + T„) < oo. □ 


Proof. Here we use the fact that for each i one has the identity deduced from the characteristic 
function of a Normal distribution, with mean 0 and variance 2 /n, evaluated at Wi = Ai — pyi. That 
is, 

/ e^ WiVi ~ TiV i ! 2 d V i = — 
v J — oo 


/ 2r i 


Now similar to the results in section 2, apply an appropriate substitution, Fubini’s theorem and the 
fact that P7 z(/3 + fyi) + nwi = N(p(/3 + fyfypi + Wifi). That is after rearranging terms it remains 
to calculate the expectation of e -JV (° n+T ”) □ 


Remark 5. Again we note that similar to the likelihood in Theorem 2.1, the likelihood incor¬ 
porating a generalized notion of leverage effects in Theorem 5.1 can be easily evaluated by classical 
numerical integration. Additionally although we have concentrated on the Barndorff-Nielsen and 
Shephard (2001a, b) models, our framework covers a large class of popular models in the literature, 
which can now be be analyzed in a likelihood framework with leverage effects. For some examples, 
see Carr, Geman, Madan, and Yor (2003).We note further that since we used an identity that does 
not depend on the distributional features of the Poisson linear functionals the results can be easily 
adapted to other processes with for instance possible additional Gaussian components. 
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